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We discuss the averaging problem in general relativity, using the form of the macroscopic gravity 
equations in the case of spherical symmetry in volume preserving coordinates. In particular, we 
calculate the form of the correlation tensor under some reasonable assumptions on the form for the 
inhomogeneous gravitational field and matter distribution. On cosmological scales, the correlation 
tensor in a Friedmann-Lemaitre-Robertson- Walker (FLRW) background is found to be of the form 
of a spatial curvature. On astrophysical scales the correlation tensor can be interpreted as the sum 
of a spatial curvature and an anisotropic fluid. We briefly discuss the physical implications of these 
results. 

[PACES: 98.80.Jk,04.50.+h] 

The gravitational field equations on large scales are obtained by averaging the Einstein equations of general rela- 
tivity (GR). The Universe is not isotropic or spatially homogeneous on local scales. An averaging of inhomogeneous 
spacetimes on large scales can lead to important effects. For example, on cosmological scales the dynamical behav- 
ior can differ from that in the spatially homogeneous and isotropic Friedmann-Lemaitre-Robertson- Walker (FLRW) 
model [1]; in particular, the expansion rate may be significantly affected. Consequently, a solution of the averaging 
problem is of considerable importance for the correct interpretation of cosmological data. It is also of importance for 
physical phenomena on astrophysical (galactic) scales. 

There are a number of theoretical approaches to the averaging problem [2-4]. In the approach of Buchert [4] a 3-1-1 
cosmological space-time splitting is employed and only scalar quantities are averaged. The perturbative approach [2] 
involves averaging the perturbed Einstein equations; however, a perturbation analysis cannot provide any information 
about an averaged geometry. On the other hand, the macroscopic gravity (MG) approach to the averaging problem 
in GR [3] gives a prescription for the correlation functions v^rhich emerge in an averaging of field equations. The MG 
approach is a fully covariant, gauge independent and exact method. We shall adopt the MG averaging approach. 
Averaging of the structure equations for the geometry of GR then leads to the structure equations for the averaged 
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(macroscopic) geometry and the definitions and the properties of the correlation tensor. The averaged Einstein 
equations can always be written in the form of the Einstein equations for the macroscopic metric tensor when the 
correlation terms are moved to the right-hand side of the averaged Einstein equations [3] . 

Spherical symmetry is of particular physical interest, and it is especially important to study the averaging problem 
within the class of spherically symmetric cosmological models. In [5] the microscopic field equations were taken and 
the averaging procedure was effected to determine the precise form of the correlation tensor in this case. In volume 
preserving coordinates (VPC), the spherically symmetric line element is given by 



ds'^ = -Bdt^ + Adr"^ 



du' 



+ 
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where the functions A and B depend on t and r . The FLRW metric in VPC is given by ( 1 ) , with A = R"^ /F'^ , B ^/R^, 
where R = R{t) and F = F{r), subject to ^ = Vl — kF^/F^ and k = —1,0 or 1. We can calculate the form of 
the Einstein tensor take averages, and obtain the appropriate form for the MG field equations and hence the 
correlation tensor C\ (for example, we have that C\ = G^f (g) — {G'^^)). In VPC, the average is then simply given by 
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which, for smooth functions with a slowly varying dependence on cosmological time, essentially reduces to a spatial 
average in terms of the averaging scale L (with L = ho/H < 1). 

The form of the correlation tensor depends on the assumed form for the inhomogeneous gravitational field and 
matter distribution. We assume that 



A{r,t) = {A{r,t)) 



1 + a„(t)L" sin (^r\ + a„(t)L" cos f^A 

n=l ^ ' n=l ^ ' 
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Bir,t) = {B{r,t)) 
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where the inhomogeneous functions A and B satisfy a set of appropriate and self-consistent conditions (for example 
{A{r,t))) = ( ^^gl''^^ /)■ Expanding in powers of L < 1, we obtain the correlation tensor up to 0{L'^) [5]: 



C" 6 = diag 



+ {AY {A), 



(5) 
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where C = r and 

e{t) = — [(ai - 36i)(ai + 61) + (ai - 36i)(ai + h)] . 

The function C can then be calculated from the contracted Bianchi identities. We note that if C" t is isotropic (i.e., 
of the form of a perfect fluid) then C = -py and (, is of the form of a spatial curvature term. Hereafter, for 
convenience we shall drop the angled brackets on averaged quantities. 

Let us first discuss averaging on cosmological scales. In the case that ^ = 0, as in the case of a FLRW background, 

the contracted Bianchi identities immediately yield C = i/A and A = 0, and £/A = £qR~^ (where £0 is a constant). 
Therefore, in this case we obtain 

6 =4i?"2^ia<?[3,l,l,l] , (6) 

and C" b is necessarily of the form of a spatial curvature term. 

The cosmological result that in the spherically symmetric case the averaged Einstein equations in an FLRW back- 
ground have the form of the Einstein equations of GR for a spatially homogeneous, isotropic macroscopic space-time 
geometry with an additional spatial curvature term, confirms the results in previous work in which we were able to 
explicitly solve the MG equations to find a correction term (correlation tensor) in the form of a spatial curvature 
[6]. This result is also (i) consistent with the work of Buchert [4], in which a spatial curvature term appears when 
averaging in a FLRW background, (ii) consistent with the results of averaging an exact Lemaitre-Tolman-Bondi (LTB) 
spherically symmetric dust model [7] , in which solutions of the LTB metric in (non-diagonal) VPC are given explicitly 
as perturbations about the spatially flat FLRW model and found to give rise to solutions which can be interpreted 
as having both spatial curvature and a constant correction term, and (iii) consistent with results in which the ef- 
fects of linear inhomogeneous perturbations perturbations on an exact spatially homogeneous and isotropic FLRW 
background [2,8] are found to give rise to correlation terms of the form of a spatial curvature term. 

Inhomogeneitics can affect the dynamics and may significantly affect the expansion rate of the spatially averaged 
"background" FLRW universe (the effect depending on the scale of the initial inhomogeneity) [2]. Therefore, a more 
conservative approach to explain the acceleration of the Universe [9] without introduction of exotic fields is to utilize 
a back-reaction effect due to inhomogeneitics of the Universe. Indeed, it has been suggested that back-reactions from 
inhomogeneitics could explain the apparently observed accelerated expansion of the universe today. This has been 
investigated by studying the effective Friedmann equation describing an inhomogeneous Universe after averaging, 
using both perturbative and qualitative analyses [8,10]. It is clear that the perturbative effect proposed always gives 
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rise to a renormalisation of the spatial curvature. It has also been argued that the effect does not simply reduce to 
spatial curvature and an acceleration can also result (although it is unlikely to be compatible with other observational 
data). 

The MG method adopted here is an exact approach in which inhomogeneities affect the dynamics on large scales 
through the correlation term. Averaging can have a very significant dynamical effect on the evolution of the Universe; 
the correction terms change the interpretation of observations so that they need to be accounted for carefully to 
determine if the models may be consistent with an accelerating Universe. Averaging may or may not explain the 
observed acceleration. However, it is clear that it cannot be neglected, and a proper analysis will not be possible 
without a comprehensive understanding of the affects of averaging. 

Let us next consider the effects of averaging on astrophysical scales (e.g., galactic scales). We assume that a galaxy 
can be approximated as spherically symmetric. In a non-FLRW background (with ^ ^ 0), the contracted Bianchi 
identities can then be integrated to obtain [5] 

^3/2^1/2. (7) 

We note that b is necessarily anisotropic (and cannot be formally equivalent to a perfect fluid). For the solution 
with 1 = (.Q = const, and 2/ = g{r)A^^, we can always write 

C b = toA-^diag [3, 1, 1, 1] - ndiag [1, 1, 0, 0] (8) 

where 11 = —{g{r)AB^/'^ — 2toA~^}. The correlation tensor C" 6 then automatically satisfies the contracted Bianchi 
identities. It can be interpreted as the sum of a perfect fluid and an anisotropic fluid (when Br 0). If both terms 
separately satisfy the contracted Bianchi identities, then the first term can be interpreted as a spatial curvature term 
and the second term can be interpreted as an anisotropic fluid with pj^ = and p|| = — Pgff. For an anisotropic fluid 
in spherically symmetric models the energy-momentum tensor is of the form diag [—^,p\\,p±,p±] , whore p|| = p+^w 
and p± = p — ^TT, and w is the anisotropic pressure. From above, we sec that if the (total) correlation tensor t is 
interpreted as an anisotropic fluid (which is comoving in VPC), it follows that 11 = — tt and p = — ^/U. Anisotropic 
fluids in spherically symmetric models have been studied in [11]. 

Although the correlation tensor C" b satisfies the contracted Bianchi identities, when interpreted as the sum 
of a spatial curvature perfect fiuid and an anisotropic fiuid through (8), the two separate fiuid do not in general 
satisfy separate conservation equations. However, the contracted Bianchi identities can be rewritten in the form of a 
conservation law for the anisotropic pressure tt, 
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in VPC where the metric is given by eqn. (1). 

Let us comment on the astrophysical applications of an anisotropic fluid. It is known that dark matter is a major 
constituent of the halos of galaxies [12]. By an analysis of observed rotation curves, under reasonable assumptions 
(e.g., that galaxies can be modelled as spherically symmetric) it has been found that the dark matter is of the form of 
an anisotropic fluid [13]. This has been taken up in [14], in which the consequences of anisotropic dark matter stresses 
are discussed in the weak field gravitational lensing (where it was noted that any attempt to model dark matter in 
galactic halos with classical fields will lead to anisotropic stresses comparable in magnitude with the energy density). 

It is of interest to further study the effects of averaging in the astrophysical context. The results of this work could 
be used to model the effects phenomenologically by including an anisotropy term (comoving in VPC), which in general 
has p = — and in the case = is of the specific form p_i = and p|| = —PqQ (where the correlation tensor is 
given by —Hdiag [1, 1, 0, 0]). The anisotropic fluid satisfies the Bianchi identities, but since it arises from an averaging 
procedure it need not satisfy any energy conditions. It may be beneficial to work in VPC, in which the metric is 
diagonal and the correlation tensor is 'comoving' (although the matter is not generally comoving). Indeed, in VPC 
the correlation tensor is given explicitly in terms of the averaged metric functions (e.g., fi = —g{r)AB^^^ — £qA~^, 
p = — 77 = g{r)AB^^^ — 2£oA~^). A disadvantage is that astrophysicists are not familiar with working in these 
coordinates. Alternatively, we could transform back to more conventional coordinates and determine the form of the 
correlation tensor; however, these coordinates may not be the most natural (e.g., the metric will not be diagonal) and 
the form of the correlation tensor (which is no longer comoving) may be quite complicated. 
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